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I. AaSTJtAC? 

Seasers to raoteljr isemsTire the lin- 
ear pclai'izatioR of crtyaae seethes bave 
been proposed for the Jfclti spectral Se- 
source Saoepler .:ntS), a xatellite sensor 
STSte* proposed to conpljjaent tbe Tbesnatic 
Mapper. At present jsstificatioc for a 
seasor oe WtS to neasare sceae pelarira- 
tioR is 1 ini ted. Tbis paper discusses a 
nodel for tbe antroat of linearly polarized 
li^t reflected by tbe sbisy leaves of 
such craps as wbeat, corn, and sorgbuiE- 
Tbe theory demonstrates that, potentially, 
•oeasureaents of tbe linearly polarized 
li^t free a crop canopy say be used as an 
additimval feature to discriminate between 
crops. Exxmiaatiaa of tbe model suggests 
that, potentially, satellite polarization 
measuresteats say be used to zonitor crop 
developeaent stage, leaf vater content, 
leaf area index, bail damage, and certain 
plant diseases. Tbe model adds to our 
understanding of the potential information 
content of scene polarization measurements 
acquired by future satellite sensor sys- 
tems such as MRS. 


II. INTROODCTION 


iaformatior, related to crop vioor and 
greartb sta'fr, inf c.-natio® need^ for var- 
ious models to predict crop grain yield*^. 
TTsere is t ccrti; .lac need to better dis- 
criminate between crops and a nee d to re- 
search and j.eveioc reiacte sensing data 
analysis tecrticues tc oredict crop grain 
yield* . 

Sensors to remotely measure the linec 
polarization of cround scenes have been 
proposed for the Kultispeetral Resource 
Sampler (KtS), a satellite sensor system 
processed to compliment the TTieisatic Map- 
per . If polarization sensors are to be 
included in the MRS, it should be estab- 
lished that polarization measurements of a 
scene provide nonredundant information ex- 
ceedinc that already found in the scene 
spectral response, now routinely measured 
by tandsat MSS with four wavelength bands 
and soon to be measured by the Thematic 
Mapper with ser'en wavelength bancs. 

At present justification for a sensor 
on MRS to measure scene p>oiarization is 
limited. Using data obtained in Lne labo- 
ratory and with an aircraft, Egan, Egan 
and Kallock, and Egan, et al. found evi- 
dence that the degree of linear oolarization 


A series of investigations has shown 
that remotely sensed Landsat satellite 
multispectral scanner (MSS) data can be 
used to accurately identify and measure 
hectarages of crops over large areas. Ex- 
perimental global wheat production fore- 
casts have been made by melding crop area 
estimates derived front Landsat data and 
estimates of crop yield made by regression 
models based on historical weather and 
yield data''^. Despite these successes, 
there are limits to the present technology; 
for example, during the Large Area Crop 
Inventory Experiment (LACIE) , there was a 
tendency for spring wheat to be confused 
with other small grains such as barley. 

And potentially in the multitemporal spec- 
tral responses of various crops there is 



Figure 1. Specular Reflection. The 
camera received specularly reflected sun- 
light from the bright areas of these 
wheat flag leaves. 
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(al (b) 

Figare 2. Polarised Light. fr» Cas- 
opT. These photos, tafcen vitb a polarizer 
oriested for transmiasiop of staziasnc spec- 
alarlr reflected light (a) aad ainlaw 
specularly reflected light (b) . denon- 
strate that the specularly reflected light 
is polarized. 

of the response cf a sceae provides addi- 
tional discriminatory iaformation with 
idsich to classify the sceaeVr*^^. Egan 
reached a potentially important conclusion 
that drying of leaves generally increases 
their depolarizing properties’*. Curran 
used a photographic measuresient technique 
to relate soil surface moisture to the 
proportion of polarized light in the scone 
response In an appendix Curran pre- 

sented data showing * possible link between 
th*: percent linear polarization of a can- 
opy and its roughness’’. 

This paper discusses a laodel for the 
ajaount of linearly polarized light re- 
flected by a plant canopy. The model is 
based upon the morphological and pheno- 
logical characteristics of the canopy and 
upon the Fresnel equations which describe 
the light reflection process at the smooth 
boundary separating two dielectrics’*. The 
theory demonstrates that the linear polar- 
ization response of the plant canopy po- 
tentially contains information to help 
discriminate between crops. The theory 
relates the response to the agronomic con- 
dition of the crop — to factors such as 
growth stage, plant vigor, and leaf area 
index (LAI) . 


III. POLARIZATION MODEL 

As may bs verified with a polarizing 
filter and camera, the shiny leaves on a 
plant (Fig. 1) provide the basis for un- 
derstanding how light is specularly re- 
flected and polarized (Fig. 2) by a 
healthy, vigorous plant stand. The shiny 
leaves of many plants, including wheat, 
corn, and sorghum vai'ieties, have a leaf 
skin or cuvicle covered by a wax layer 
which spc'^nlarly reflects light in accor- 
dance with the Fresnel equations. Because 


the light is specularly reflected, the 
equations show tiaat it :s polarized for 
all bet two angles of iacideace (0 aad 9S 
degrees). Brrece and Bclaes fooad that 
the bidirectional light scattering charac- 
teristics of wheat, com. and soybean 
leaves istermeciate to t3ae scattering 
diaracteristics of diffuse and specular 
reflectors, indicating that the specular 
portion is an important part of the total 
leaf response”. 

Proa: the laws of optics a leaf seg- 
ment with a tendency toward specularity 
portends that a significant portion of the 
surface area is flat and similarly orient- 
ed. Indeed, there is electron microgra{^ 
evidence that the wax deposits on glossy 
leaves may form saocth films on the cuti- 
cle or platelets which lie flat on the sur- 
face ’ . Electron microgra^s of a wheat 
leaf and a com leaf both reveal irregular 
acicular wax structures distributed on a 
flat wax surface much like tree stumps on 
a flat, clear cut area or a child's jacks 
scattered ok a table ’*•**. 

A. THEXJRT 

The mathematical model for polariza- 
tion of light from a idieat canopy is de- 
veloped in two parts. First, the micro- 
scale situation, the pclarizatlon response 
and orientation requirements for a small 
specularly reflecting ares ia on a leaf 
(Fig. 3), is analyzed. Second, the micro- 
scale results are extended t. the macro- 
scale level as measured by a field spec- 
trcHseter or satellite sensor (Fig. 4). 

The assumptions are: 

1. There exist on the wax surface of 
each shiny leaf small flat areas, £a, 
which specularly reflect light. 

2. The wax layer is essentially clear 
and absorbs little light. This means 
that for the wax layer the complex index 
of refraction can be adequately approxi- 
mated by its real component, a reasonable 
supposition for the visible spectral re- 
gion where any light energy absorbed by 
the wax layer is then unavailable to the 
chloroplasts to promote photosynthesis. 
Limited evidence supports tliis assumptiorf®. 

3. Specular light reflection occurs 
principally at air-wax boundary. Com- 
paratively negligible amounts of light are 
reflected specularly to an observer and 
from the boundaries between epidermal cell 
walls, cell membranes, and the various 
cutible layers. These boundaries have 
comparable indices of refraction and often 
appear rough in electron micrographs . 

4. The magnitude of polarized light frexn 
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taoise •oil. Cor exotic) otter 
lit len«* i« iAsi^oificaat. 


Wgoro Scale tespoaoo . Ob the micro 
•calc'Teeal, Pis- 3, ao^i^t i* »pecelar- 
ly reflected by <>■• of the aaall leaf 
areas to aa ofaserver oely if tte area is 
iweperly orieeted. the eoreal to the area 
east be la tte plane and bisect the mqle 
Coned by tno sectors, tte illninatioa 
s e ct or (E) directed Cron the leaf ar^ to 
tte san and tte afasersatian sector (Vi di- 
rected Cron tte leaf area to tte observer, 
m otter eords tte area £a east be orient- 
ed son that tte anyle of incidence, y, 
eqeal^s Ite angle of reflectance and such 
that tte sectors E.^Sa (the imit sector 
nomal to £a) , and V are ooplanar. Only 
then vill specularly reflected li^t Cron 
Aa rean an obserser. these con<*l'' ons 
fore tte keystone of tbe polari: ^ ->a 

nodel and are satisfied by tte e .ations. 


®-"a 

Y*COS * r *cos* 

n 



=0.5cotf> 


the solar irradiance incident on one 
snail area, Aa. is a function of the angle 
of incidenre. If the area Aa specularly 
reflects light to an observer, then the 
angle of incidence is unicmely determined 



X toward 
sun 


Figure 3. Coordinate System. A 
small leaf area Aa specularly^ref lects 
sunlight toward an o^sejver, if and 
only if the vectors E, n^ and V are co- 
planar and the angles of incidence (y) 
and reflectance are equal. 




h 


Figure A. Canopy Response. A sensor 
measures the canopy response over a solid 
angle Aw. 

by tbe angles ( 65 . 9^'^'* ^ discussed 
above. For such an area the radiant flux 
Incident on Aa is Psl*[ cos y %fhere y = 

prteability of 
finding in a small volume an area Aa il- 
luminated directly by the sun as opposed 
to being shaded by intervening foliage. 

The probability of finding in a small vol- 
ume an observable area Aa is symbolized by 
Pv In determining Ps and P^. area Aa is 
assumed to be either illuminated or shaded, 
observable or not observable; adundaral 
effects are not considered**. 

The probabilities, Pg and Py, are 
functions both of the (x,y,z) location of 
the leaf in the canopy and of directions 
of illumination and observations, respec- 
tively. The probabilities will be unity 
when only leaves are illuminated and ob- 
served. For example, Pg and will ap- 
proach unity for the topmost leaves of a 
dense, preheaded wheat canopy if the ag- 
gregation of these leaves forms a layer 
one leaf thick at the extreme top of the 
canopy essentially impenetrable to direct 
illumination. The probabilities will be 
less than unity for more typical canopies 
with some soil and/or non-leaf foliage 
illuminated and observed and some leaves 
not illuminated and/or not observed. 

Even though the incident sunlight is 
not polarized, each small area, Aa, polar- 
izes the specular portion of the reflected 
sunlight provided the angle of incidence 
is neither 0 or 90 degrees. If Aa is 
smooth, the magnitude of the light that 
is specularly reflected and polarized by 
Aa is described mathematically by the 
Fresnel equations and Stokes vector and 
depends only upon the angle Of incidence, 
Y, and the index of refraction of both 


100 


1980 Alachine Processing of Rerrx>lely Sensed Data Syirposium 


ORIGINAL V:\G 
OF POOR (M ’ A 


tiK epicoticale wax lay«r wd air**. 

Vtor a lapotfa svrface the portions of 
tbe specular reflectance with electric vec- 
tor perpendicular. Pspl^Y). eod parallel, 
^?p! i *** plane of incidence are 
given 1^ the Fresnel equations**. Ihe 
first oowponent of the Stokes vector, Sx# 
provides the waqnitade of light reflected 
hf a surface**. The tern $x is tlte son of 
the specularly (Sg) and diffusely (Sq) re- 
flected light fron the surface. The sec- 
ond ooaiponent, Sq, is the portion of Sx 
linearly polarized by tile surface. 


®I * 


(2a) 

Ss = 


(2b) 


“^spi-^’spll’/^.O 

(2c) 


[dinensionless] 


Most often Aa is not a perfectly 
snooth surface but instead supports snail 
acicttlar structures. These structures 
diffusely scatter light which would other- 
wise be specularly reflected. To account 
for the reduced anount of light, specular- 
ly reflected by an area &a idiich is not 
perfectly snooth, the Fresnel equations 
is modified by a factor, K. In general K 
[dimensionless] is a function of the 
wgles (6s,9y,4v)« wavelength, and side 
(top or bottom) lateral position and di- 
rection on the leaf. However, evidence 
suggests that for any one leaf the wax 
acicular structures in most cases are ho- 
mogeneously and isotropically dispersed 
across the leaf surface”. It is assumed 
here that K(9s,ev«<^v'^) is identical for 
all leaves and is not a function of lat- 
eral position and direction on the leaf 
surface. The value of K varies between 
zero and one. 

Define a probability density function 
fa(6,i(i) for leaves such that the probabil- 
ity that any one of the leaf areas Aa is 
oriented within a solid angle Aug about 
(9a»<(’a) is A(Oafa(ea,^a> • 

The units of fa=fa(x,y,z,6,^) are [sr~*]. 
Because the area Aa must be correctly ori- 
ented to reflect light to an observer, the 
Jacobian provides*’ 

A(i)g = Aui^/4cosY = AyCOS*9^4h*cosY (3) 


area is tbe quadnqile product of (1) the 
radiant flux incident on the area Aa, (2) 
the specular reflectance of Aa, (3) the 
probability that Aa is correctly oriented 
to specul;^ly reflect li^t in a solid 
angle Aw^ .^boot directi<» (€v«9v) , *od (4) 
the probability that Aa is observed. 


(1) (2) (3) (4) 

"^♦as * 

For a volume Vj, = ||«**aS <5) 

all leaf area in Vj 

If the volume Vj is sufficiently sasall, 
then the probabilities Ps and Pv are essen- 
tially constants everywhere in Vj and will 
be denoted Pgj and P^j- Letting Aj be the 
leaf area in volume Vj, f^^j be the proba- 
bility density function of leaf area ori- 
entation in Vj, and using eq.i3) to sub- 
stitute for Au^, eg. (5) becomes 

*Vj=Psj !E|AjKSsAvCOs*6vfaj (9a,4a) Pvj/<h* (6) 


Sunming the specular flux contributions of 
each volume Vj in the field of view of the 
sensor to find the specular portion of the 
flux measured, eq. 7, and the linearly po- 
larized portion of *s, eq. 
all Vj in 
field of view 


1 

♦S=|E| KSg, 

^ ^ajPsj’^v^^j 

Aycos’e v'4h* 




source 

dependent 

term 

canopy 

dependent 

tepis 

sensor 

dependent 

terms 


all Vj in 
field of view 


4Q=|Ej KSq Z^fajPgjPvjAj Aycos’e v/4h* 


(7) 


( 8 ) 


The percent linear polarization is propor- 
tional to the linearly polarized flux di- 
vided by the total flux, the sum of the 
diffuse and specular fluxes. 


100%$q/($d + *s) 


Macro Scale Response . On the macro 
scale level. Fig. 4, the radiant flux due 
to specularly reflected sunlight received 
from leaves in the field of view of a sen- 
sor is found by summing the flux contrib- 
uted by each leaf area Aa in each volume 
Vj in the field of view. For a Aa the 
radiant flux specularly reflected into a 
solid angle Aw^ by a randomly selected 


To illustrate the properties of the polar- 
ization model, the response of several 
canopies will be examined. 

Example A: Sparse Wheat Canopy . If 

the properties of the canopy are constants 
for those layers containing leaves, that 
is, if Psj“Psi“Ps' Pv-i"Pvi=Pv» Aj=Ai=A, 
and faj=fai=fa' ^Vj,V£E field of view 
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Figure 5. Preheaded Wheat Canopy 
Polarizaticn Response. Prior to heading 
the response is zero at the anti-solar 
point, the "hot spot," and increases with 
increasing ze.iith view angle. 



amlch Mfl* 

AaJaatk vlov mmmf fcM mm ■ Asiawtk rim trm4 mm 

Figure 6 . Headed Wheat Canopy Polar- 
ization Response. After heading the re- 
sponse remains zero at the anti-solar 
point, is maximum at intermediate zenith 
view angles, and approaches zero for near- 
horizontal view directions where heads and 
stems obstruct view of polarizing flag 
leaves. 


and 

•s=iE|Ps(LM)KSsfaPvAvAoi/2cos*ev (10a) 

^ m » 

source canopy sensor 



♦0="!E|Ps(IAI)KSQfaPv AwAv/2cos9v UOb) 

where LAI is leaf area index. Ttie linear- 
ly polarized portion of the total scene 
radiance is*\ 

Lq=*q/A^=^| PsUVIKSQfaPv/2cos6v (11) 

Example B: Preheaded, Dense Wheat 

Canopy . If the probability PsjPvj'l 
Vj in the topmost layer of the canopy and 
PgjPvj^O for Vj in all lower layers, then 
the ®re proportional to the leaf 

area index only of the topmost layer. A 
winter %fheat canopy measured just prior 
to heading might have the following char- 
acteristics: LAI =2.0 for top layer con- 

taining flag leaves with a wax layer index 
of refraction = 1.5, fa = uniform = 0.0796 
sr~*, K = 0.9, and Pg = Py = 1.0 for top 
layer and P- = Py = 0.0 for all lower lay- 
ers. The linearly polarized portion of 
the canopy radiance (eq. 11) and the lin- 
early polarized flux measured by a sensor 
over such a canopy are shown in Pig. 5. 

The calculations are for a sensor with a 
field of view, Au, of 15“ = 0.216sr, en- 
trance optics of area 0.002m^, and spec- 
tral band of 0.6-0.7um (red '-'avelengths) . 
The solar insolation is assumed to be 
165.3W/m^ in the 0.6-0.7um spectral band^^. 

Example C: Headed, Dense Wheat Can - 

opy . The pol::rization response of a wheat 
canopy is expected to change significantly 
during the heading growth stage (Fig. 2) . 
This is because the probabilities Pg and 
Py of a headed wheat canopy, unlike those 
of a preheaded canopy, are pronounced 
functions of sun angle and view angle. 

The product Ps^v *’'®y estimated for a 
hypothetical headed canopy with LAI =2.0 
by applying linear regression techniques 
to data for a canopy with LAI =1.0 and 
scaling by a factor of 2,0^*. 

(12) 

PsPy=.232exp(-1.17 ( (l/cos6y) + (l/cosOg) ) ) 


s.t.A?«0, then 

all Vj in 
field of view 

^ “sj^vj^j 
j=l ^ J 

= PgPy[total two sided leaf area in POV] 
= PsPy[2LAl][h*Aut/cos’eu)y] 


Equation 12, derived assuming Pg and Py to 
be independent, provides erroneous esti- 
mates of PsPv angles near the canopy 
hot spot direction where the probability 
Pdeaf observation I leaf illumination) 
approaches unity. Fig. 6 shows the linear 
polarization response for a source-canopy- 
sensor with the parameters of Example B 
except PgPy given by equation (12) . 
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Example Pi Percent ^lari««tiy 
RespoSe ot Cay>py . xt the nagnituSe of 
the specular flux is either saall or large 
cocapared to the diffuse flux, then the 
percent polar ixation (eq. 9) reduces to 

{ 100%*q/*d.*o»*s (13a) 

100%*q/«s**D<"*S 

Both 9 q and are proportional to the 
saw agronoaiic factors. When 9 q<<*s 
( idiich Bight occasionally be true in the 
blue and red spectral regions), then there 
is no agronomic inforaiation in the term 
percent polarization because the agronomic 
factors in the numerator and denominator 
of eq. 13b cancel. 


IV. DISCUSSION 

The model shows the magnitude of the 
polarization response of a plant stand is 
related to the solar insolation and the 
characteristics of the canopy and the sen- 
sor. The response depends on the optical 
and geometric properties of the portion of 
the canopy in the instrument field of view. 

The calculatloi.s show that a sensor 
would measure zero linearly polarized 
light at the anti-solar point, the canop'^’ 
"hot-spot," where Sq and the angle of in- 
cidence of the sunlight are both zero. 

The sensor would measure the maximum 
amount of polarized light in the solar 
azimuth direction (Fig. 5), provided the 
small areas &a are randomly oriented in 
azimuth and zenith directions. Otherwise 
the direction of maximum polarization may 
be shifted, as might occur when a strong 
wind preferentially orients the flag 
leaves of wheat downwind. 

The theory shows that when the ap- 
proximation (eq. 13b) is valid — when the 
specular flux is much, much greater than 
the diffuse flux — the percent linear po- 
larization is not directly related to the 
canopy agronomic properties. This might 
occur in the chlorophyll absorption re- 
gion in the red portion of the spectrum 
viewing at large zenith angles toward the 
sun azimuth angle. However, even though 
the percent linear polarization in certain 
circumstances may contain limited infor- 
mation related to canopy agro.nomic fac- 
tors, the magnitude (eq. 8) of the linear- 
ly polarized flux is always directly re- 
lated to the canopy agronomic properties . 
Thus, the model provides a theoretical 
basis for the same, but empirically based 
result noticed by Egan’ . 

The agronomic variables LAI, fg, K, 

Sg, Pg, and Py in the equations (eqs. 8 


and 10b) are functions of one or more en- 
vironmental and/or physiological variables. 
The leaf area index (LAI) is a function of 
many variables including species and cul- 
tivar, weather, and growth stage. The 
prc^ability density function fa(9a<9a) 
the orientation of the leaf areas Aa is a 
function of wind strength and direction, 
catastrophes (such as hail damage) , crop 
vigor (moisture stress causes corn leaves 
to roll , cotton leaves to droop) , and 
growtii stage (the shapes of healthy green 
leaves and senescent leaves are not the 
same) . The optical properties of the 
leaves, K, Sg, and Sj, are functions of 
species and almost certainly cultivar, di- 
sease (plant pathogens often alter or de- 
stroy the wax layer) , pubescence (the 
hairs scatter light which would otherwise 
be specularly reflected) , material on the 
leaves (dust, pollen, water droplets), and 
wavelength (the Fresnel equations are func- 
tions of the cuticle wax index of refrac- 
tion, which changes in regular fashion 
with wavelength) . The terms Sg and Sg are 
functions of the angle of incidence (view 
and illumination direction) of the sun- 
light. The geometrical properties of the 
canopy, Pg and Py, are functions of the 
angles of illumination and observation 
(lower leaves have lower probabilities of 
being illuminated and observed at large 
sun and view zenith angles) and growth 
stage (wheat heads partially block illumi- 
nation and observation of flag leaves; the 
projected ar>‘>a of leaves changes with se- 
nescence) . Se 'eral of these functional 
relationships vill be discussed further 
in the folxovrng paragraphs. 

Light polarized by a moist or wet 
soil surface is a part of the canopy po- 
larization response not considered in the 
theory. Visual observations suggest that 
except for wet soil surfaces, the amount 
of light polarized by the soil is insig- 
nificant compared to the amount of light 
polarized by foliage. Neglecting emergent 
and sparse canopies, the soil generally 
has a very low probability of being both 
illuminated and observed. Thus, the the- 
ory presented predicts the polarization 
response of canopies on dry soils and/or 
with sufficient foliage to obscure the 
soil . 

Sunlight tends to be specularly re- 
flected and polarized by leaves in the 
upper portion of plant stands. The proba- 
bility that a leaf is both visible and 
directly illuminated by sunlignt is often 
a pronounced function of the (x,y,z) lo- 
cation of the leaf in the stand and of 
the view and illumination directions. 

The probability tends to decrease rapidly 
with increasing depth into the canopy. 

This means that leaves in the lower portions 
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of a plant stand vill little affect the 
canopy polarization response. Thus, lower 
leaf senescence or a disease condition lo- 
calized to the lower leaves may not be de- 
tectable using polarization measurements. 

The probability density function, 
fa^^st^a)' orientation of the leaf 

areas, Aa, can be calculated (from eg. 10b) 
for tte population of observable, specu- 
larly reflecting leaves using polarization 
measurements sampled from the hemisphere 
of all possible canopy view directions. 
Such a density function is needed as input 
data to certain canopy radiation models 
which are used to ex^unine the utility of 
and information in canopy reflectance 
measurements*^. To obtain the polarization 
data needed to calculate f^, Horvath pro- 
posed a field apparatus consisting of a 
linearly polarized light source co-located 
with a sensor with a polarization analy- 
zer**. Due to inherent randomness of 
the leaf strut; ere, the diffuse portion 
of the reflected light will '.end not to be 
polarized'*. The specular portion of the 
reflected light will be polarized, not 
because of reflection at the leaf cuticle 
wax layer (0® angle of incidence] but in- 
stead because the light source is polar- 
ized. To compute f^ the canopy polariza- 
tion response is measured with the analy- 
zer oriented in two directions, parallel 
and perpendicular to the polarization di- 
rection of the light source beam. It is 
easily shown that the specularly reflected 
flux (eq. 10a) is the difference of these 
two measurements; hence, in a particular 
direction 

fa = |e|LAI KSsP3pvawAy/2$sCOS6v 

where Sg = ( (n-l)/(n+l) ) * for normal inci- 
dence, and Pg=Py. A practical limitation 
to the approach exists. In general values 
for leaf area index fLAl) , the factor K, 
and the index of refraction (n) , are not 
known; however, properly normalizing to 
unity the integral with (Oy/ijiv) of the ini- 
tial estimate of fg obviates the need to 
know these terms. But more importantly 
the need to know or estimate P^ and P^, 
usually unknown functions of view direc- 
tions, cannot be circumvented if fa is to 
be calculated. The term PgPyfa includes 
all the canopy dependent variation due to 
view angle (assuming the factor K is a 
constant) and is always calculable. 

Depending on its direction and 
strength, the wind is capable of reorient- 
ing the leaves of a canopy and thereby 
changing the probability density function 
of leaf area, fg. The resultant varia- 
tions from day to day in the polarization 
response of a field will tend to compli- 
cate interpretation of polarization data 


because these variations represent noise 
(unless, of course, the probability densi- 
ty function of the orientation of the 
leaves is to be calculated for each day) . 
The size of these day to day variations, 
if sufficiently large, might preclude a 
naive analysis %diich neglects wind effects; 
conceivably, wind induced variations in fa 
might render an agronomic interpretation 
of polarization data in^ssible. There 
remains the hope, however, that for data 
taken at one time the wind will affect 
similarly all the fields in a region con- 
taining a particular crop species and cul- 
tivar. For such data comparisons between 
fields of a specific crop variety might 
remain valid. Therefore, an important 
question which should be addressed empiri- 
cally is the following: How uniformly 

does wind affect polarization data acquired 
over a region? 

Visual observations suggest that blue 
skylight incident on the canopy affects 
minimally the magnitude of the canopy po- 
larization response. By the same process 
discussed in the theory for sunlight, the 
shiny leaves of a canopy polarize the 
specularly reflected skylight, a spectrally 
varying light source already polarized 
according to observer view direction. The 
magnitude of the skylight and its effect 
on the canopy polarization response is 
greatest in the blue spectral region and 
decreases into the near infrared. Atmos- 
pheric haze, which decreases the solar 
insolation on leaves, noticeably decreases 
the specular and polarized light from 
leaves . 

The efficacy and feasibility of a 
satellite sensor measuring the linear po- 
larization of a sc-^ne through the earth's 
atmosphere has not been considered. Sen- 
sor design must consider the path radiance 
of the atmosphere, a source of linearly 
polarized light potentially capable of 
altering or masking the amount of linearly 
polarized light received from the scene. 
Even if the polarized portion of the path 
radiance is excluded, analysis of field 
spectral radiometer data suggests that for 
a satellite polarization sensor the signal 
flux must be increased (the spectral/spa- 
tial resolution of a Landsat-type satel- 
lite must be degraded probably by a factor 
greater than 10:1) and/or noise power de- 
creased to obtain a signal to noise ratio 
approximately equivalent to that of Land- 
sat in the red spectral channel. Including 
atmospheric effects in the analysis would 
potentially indicate a practical value of 
the spatial/spectral resolution of a sat- 
ellite sensor. 

The information in canopy polarization 
data, when obtained from satellite sensors, 
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(a) (b) ic) 


Figure 7. Polarization Response with Heading. The three photographs, taken at one 
week intervals before and after headirg, demonstrate that the amount of specularly re- 
flected- and therefore polarized-sunlight decreases significantly with heading. 


probably will be used in conjunction with 
other remotely sensed data and will be 
extracted by analysis of frequent, synop- 
tic data sets, by using the temporal and 
spatial information to make relative com- 
parisons between the fields in the data 
for one date and between the dates for 
one field. One polarization measurement 
of one field for one date probably will 
have little value unless it is compared 
to polarization data for that field and 
other fields for that date and other dates. 
This is because it is unrealistic from the 
model to expect that canopy polarization 
data will be calibrated in an absolute 
sense to discriminate a particular crop 
or to correlate uniquely to a particular 
agronomic variable. Frequent, synoptic 
polarization data from a satellite sensor 
potentially aid in assessing crop vigor 
and growth stage and in determining areas 
of hail damage and pestilence, all poten- 
tially possible from comparisons between 
field and across dates. Perhaps daily 
satellite coverage is feasible using a low 
spatial and spectral resolution sensor in 
a geosynchronous orbit. 

The canopy polarization response de- 
scribed by the model is a function of 
wavelength only because the index of re- 
fraction of the cuticle wax layer is a 
function of wavelength. From the physics 
of the optical propertiec of materials it 
is expected that the index of refraction 
of the wax layer will gradually and mono- 
tonically increase with decreasing visible 
wavelength, displaying no perturbations 
or "fine structure" with wavelength''*. 
However, the model indicates the percent 
linear polarization of a healthy green 
canopy will be large in the blue and red 
spectral regions, small in the green, and 
even smaller in the near infrared region 
away from any absorption bands. This is 
because the total canopy flux, the normal- 


ization factor used when computing percent 
polarization, exhibits a green vegetation 
response^’. 

From the model there appears little 
need to measure the canopy spectral polar- 
ization response with high wavelength re- 
solution in the visible spectral region; 
a polarization sensor covering the entire 
visible region or a large portion of it 
might suffice. Conversely, in the infra- 
red spectral region the cuticle wax layer 
may absorb in narrow spectral regions de- 
fined by the structural properties of the 
constituent waxes of the layer, by the 
resonant frequencies of the translational 
and rotational vibration modes of mole- 
cules of the layer''*. If absorption bands 
exist, high resolution spectral polariza- 
tion data may possibly provide information 
concerning the properties of the cuticle 
wax layer, properties relatable to crop 
species and light regime'®. 

The linear polarization model may be 
extended to include the elliptical polar- 
ization response of the canopy. Evidence 
exists that the cuticle wax of some species 
is birefringent and therefore potentially 
capable of elliptically polarizing specu- 
larly reflected light'. However, this 
evidence does not suggest that the struc- 
ture of the wax layer is sufficiently 
organized to elliptically polarize light 
specularly reflected from a significant 
amount of leaf area. Egan argues that 
the amount of elliptically polarized light 
from a vegetation scene should be negli- 
gible because of the inherent randomness 
of the properties of the vegetation’ 

The connection between leaf polariza- 
tion measurements and leaf moisture con- 
tent, noted by Egan, is supported by a 
morphological model for the structural 
changes which occur in a leaf undergoing 
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dehydration". When the leaf water content 
and leaf thickness decrease and the leaf 
cells dehydrate and collapse, the surface 
roughness of the wax layer increases. The 
leaf surface appears rough in thin sections 
of senescent leaves and of leaves under 
moisture stress**. As the surface rough- 
ness increases, the specular portion of 
the light reflected by the leaf decreases 
because there are fewer areas Aa which are 
similarly oriented. The amount of linear- 
ly polarized light reflected by the canopy 
decreases in company with the decrease of 
specularly reflected light. These argu- 
ments suggest that the canopy polarization 
response should decrease with decreasing 
leaf water content in the canopy and 
therefore serve as an indicator of canopy 
moisture stress. Visual evidence supports 
this hypothesis. Leaves under moisture 
stress often appear less shiny than fully 
turgid leaves. Dry, senescent leaves 
often have a matte surface finish. 

Detection of the date of heading of 
a wheat canopy (Fig. 2), information which 
is needed for use with phenologically 
based models to predict the ultimate grain 
yield of the crop, might be feasible using 
satellite polarization measurements (Pig. 

7) . The eventual weight of grain produced 
by each wheat plant is largely determined 
by the condition of the flag leaf, its 
size and vigor, and by the weather regime 
endured by the plant following heading 
when the grain head begins to fill**. Know- 
ledge of the date of heading permits a 
better estimate of the post-heading weather 
for the crop. Prior to heading the top- 
most foliage on the wheat plant is the 
flag leaf, easily the most visible and 
Illuminated canopy component (Fig. 2) . 
Following heading, wheat heads are the 
topmost foliage and partially obscure the 
flag leaves to both sunlight and observa- 
tion, changing the values of both Ps and 
Pv Figures 5 and 6 show that the magni- 
tude of the polarized light, which depends 
directly upon the specular reflections 
from flag leaves, will decrease by a fac- 
tor of 60 for ©5=30 and 6y=0 for the two 
hypothetical canopies during heading as 
the leaves are increasingly obscured to 
both illumination and observation. The 
obscuration of the flag leaves is enhanced 
at off nadir observation angles directed 
toward the solar azimuth (Fig. 6) . Poten- 
tially both the condition of the flag 
leaves and the date of heading of a crop 
might be monitored using polarization 
measurements obtained from a satellite 
sensor with both on and off nadir viewing 
capability. Such a view capability has 
been proposed for the MRS sensor. 

Applicability of the polarization 
model should extend to many species because 


shiny leaves trtiich specularly reflect sun- 
light are ubiquitous, unconfined by geog- 
raphy or climate. Other plants besides 
wheat, sorghum, and corn with specularly 
reflecting leaves include coffee, Sudan 
grass, banana, orange, sugarcane, and many 
forest species. Schieferstein and Loomis 
found epicuticular wax deposits on about 
half of the plant species they tested'*. 
However, the mere presence of a cuticle 
wax layer does not guarantee that a leaf 
will specularly reflect and polarize a 
significant portion of the incident light; 
the leaf must also appear shiny. Fibrillar 
light scattering significantly diminishes 
the polarization response of pubescent 
soybean leaves. And the surface of the 
wax layer of some species is insufficient- 
ly smooth to specularly reflect light". 


V. CONCLUSIONS 

This paper discusses a model for the 
amount of linearly polarized light reflect- 
ed by the shiny leaves of such crops as 
wheat, corn, and sorghum, each a grain of 
major economic importance to the world. 

The model is based upon the morphological 
and phenological characteristics of the 
canopy and upon the Fresnel equations 
which describe the light reflection pro- 
cess at the smooth boundary separating two 
dielectrics. 

The theory demonstrates that, poten- 
tially, measurements of the linearly po- 
larized light from a crop canopy may be 
used as an additional feature to discrim- 
inate between crops such as wheat and 
barley, two crops so spectrally similar 
that they are misclassif ied with unaccep- 
table frequency. Examination of the model 
suggests that, potentially, satellite po- 
larization measurements may be used to 
monitor crop development stage, leaf water 
content, leaf area index, hail damage, and 
certain plant diseases. Such information 
is needed for use with models which pre- 
dict crop grain yield. 

The model adds to our understanding 
of the potential information content of 
scene polarization measurements. The in- 
formation content of these measurements 
has not been extensively investigated and 
needs to be understood to evaluate the 
potential usefulness of the proposed po- 
larization sensor for the satellite borne 
Multispectral Resource Sampler. The effi- 
cacy of a satellite sensor measuring the 
linear polarization of a scene through the 
atmosphere remains to be determined. 
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